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Inclusive multi-jet production is studied using the ATLAS detector for proton-proton collisions with a center- 
of-mass energy of 7 TeV. The data sample corresponds to an integrated luminosity of 2.43 pb~^, using the first 
proton-proton data collected by the ATLAS detector in 2010. Results on multi-jet cross sections are presented 
and compared to both leading-order plus parton-shower Monte Carlo predictions and next-to-leading-order QCD 
calculations. 



^ 1. Introduction 

^ . 

O ■ 

CN ' . . n 

This document represents a summary of the analysis and results presented in 

Q-i' 
^ '. 

m , At hadron colliders, events containing multiple jets in the final state are prolific and provide a fertile testing 
ground for the theory of the strong interaction, quantum chromodynamics (QCD). 

■ In this paper, a first study is performed of multi-jet events produced from proton-proton collisions at 7 TeV 
' center-of-mass energy, using the ATLAS detector at the Large Hadron Collider (LHC) at CERN. The data 

^ . sample used for this analysis was collected between April 10 and August 30 of 2010 and represents a total 

■ integrated luminosity of 2.43 pb"-'^. 

^ ^ The two primary motivations for the multi-jet study in this paper are to evaluate how robust leading-order 
^ (LO) QCD calculations are in representing the high jet multiplicity events, relevant as backgrounds in high- 
ly . energy searches, and to test next-to-leading-order perturbative QCD (NLO pQCD) calculations. 

> 

cn , 2. The ATLAS Detector 

00 

• ' The ATLAS experiment consists of an approximately 45-meter long, 25-meter diameter cylindrically shaped 
I detector centered on the proton-proton interaction point. A detailed description of the ATLAS experiment can 

■ be found elsewhere [3]. 

High-energy particles produced in collisions initially pass through an inner tracking system embedded in a 
2 T solenoidal magnetic field. Just outside the inner tracker package is a system of liquid argon and scintillating 

. , , , tile calorimeters used for the measurement of particle energies. A liquid- argon/lead electromagnetic calorimeter 

covers the pseudorapidity^ range of |r/| < 3.2. This calorimeter is complemented by hadronic calorimeters, built 
H using scintillating tiles and iron for \r]\ < 1.7 and liquid argon and copper in the end-cap (1.5 < \r]\ < 3.2). 
. . . 1 Forward calorimeters extend the coverage for ATLAS to \r]\ = 4.9. The calorimeters are the primary detectors 
used to reconstruct the jet energy in this analysis. 

The ATLAS trigger system employs three trigger levels, of which only the hardware-based first level trigger is 
used in this analysis. Events are selected using the calorimeter-based jet trigger. The first level jet trigger Q uses 
coarse detector information to identify areas in the calorimeter where energy deposits above a certain threshold 
occur. A simplified jet finding algorithm based on a sliding window of size A0 x At; = 0.8 x 0.8 is used to 
identify these areas. This algorithm uses coarse calorimeter towers with a granularity of Ac/) x Arj ~ 0.2 x 0.2 
as inputs. 



ATLAS uses a right-handed coordinate system with its origin at the nominal interaction point (IP) in the centre of the 
detector and the z-axis along the beam pipe. The x-axis points from the IP to the centre of the LHC ring, and the y axis points 
upward. Cylindrical coordinates (r, 0) are used in the transverse plane, (p being the azimuthal angle around the beam pipe. The 
pseudorapidity is defined in terms of the polar angle 6 as rj = — lntan(0/2). The rapidity is defined as j/ = 0.5 X ln[(_E+pz)/(E— pz)], 
where E denotes the energy and pz is the component of the momentum along the beam direction. For massless objects, the rapidity 
and pseudorapidity are equivalent. 
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3. Cross Section Definitions and Kinematics 

The anti-fcf algorithm 0, Q with full four-momentum recombination is used to identify jets. For high mul- 
tiplicity studies, namely up to six jets in an event, the resolution parameter in the jet reconstruction is fixed 
to i? = 0.4 to contend with the limited phase space and to reduce the impact of the underlying event 6] in 
the jet energy determination. For testing NLO pQCD calculations, where the study focuses on three-jet events, 
a resolution parameter of i? = 0.6 is preferred, since a larger value of R is less sensitive to theoretical scale 
uncertainties. 

Jet measurements are corrected for all experimental effects and refer to the particle-level final state. At the 
particle level, jets are built using all final-state particles with a proper lifetime longer than 10 ps, including muons 
and neutrinos from hadronic decays. These corrections are described in Section [B) The NLO pQCD calculation 
is not interfaced to a Monte Carlo (MC) simulation with hadronization and other non-perturbative effects. The 
correction for non-perturbative effects applied to the NLO pQCD calculation is described in Section H) 

Cross sections are calculated in bins of inclusive jet multiplicity, meaning that an event is recorded in a jet 
multiplicity bin if it contains a number of jets that is equal to or greater than that multiplicity. 

Inclusive multiplicity bins are used because they are stable in the pQCD fixed-order calculation, unlike 
exclusive bins. Only jets with > 60 GeV and \y\ < 2.8 are counted in the analysis. These cuts are chosen 
to ensure that the jets are reconstructed with high efficiency. The leading jet is further required to have 
PT > 80 GeV to stabilize the NLO pQCD calculations Q- 



4. Theoretical Predictions 

Measurements are compared to pQCD calculations at leading order and next-to-leading order. For the leading- 
order analysis, ALPGEN 8] is used to generate events with up to six partons in the final state using the leading- 
order set of proton pdfs CTEQ6L1 ^. ALPGEN is interfaced to PYTHIA 6.421 [13, [HI and, alternatively, to 
HERWIG/ JIMMY [12Hl5j to sum leading logarithms to all orders in the parton-shower approximation and to 
include non-perturbative effects such as hadronization and the underlying event. The ATLAS generator tunes 
from 2009 (MC09') ^ and from 2010 (AUETl) [l^ are used. Additional tunes are investigated to assess 
the impact of the underlying event and parton-shower tuning. With comparable underlying event tunes and 
ALPGEN parameters, the comparison between ALPGEN-I-PYTHIA an ALPGEN-f HERWIG/ JIMMY uncovers 
differences that may arise from different parton-shower implementations and hadronization models. 

SHERPA [18,] with its default parameters and renormalization scale scheme from version 1.2.3 is also used 
to generate events with up to six partons in the final state. This provides an independent matrix element 
calculation with a different matrix element and parton shower matching scheme. 

The PYTHIA 6.421 [13, [HI event generator is included to study the limitations in a 2 2 calculation. The 
MRST 2007 modified leading order pdfs interfaced with the AMBTl[i3| generator tune is used in the 

sample generation. 

The generated particles are passed through a full simulation of the ATLAS detector and trigger [2l| based on 
GEANT4 ^2^] to account for detector effects. Additional proton-proton collisions are added to the hard scatter 
in the simulation process to reproduce realistic LHC running conditions. Events and jets are selected using the 
same criteria in data and Monte Carlo simulations. 

For the NLO pQCD study, the calculation implemented in NLOJet-|--|- [2^ is used. The renormalization and 
factorization scales were varied independently by a factor of two in order to estimate the impact of higher order 
terms not included in the calculation. An additional requirement that the ratio of the renormalization and 



factorization scales did not differ by more than a factor of two was imposed. Two NLO pdf sets CTEQ 6.6 [24 1 
and MSTW 2008 nlo [13] were used for calculating the central values. The 90% confidence-limit error sets are 
used in the evaluation of the pdf uncertainties. The uncertainty in the measurement due to the uncertainty in 
the value of as is calculated by varying the value of as by ±0.002 for each pdf set. 

The NLOJetH — h program is a matrix element calculation, and therefore it lacks a parton-shower interface and 
does not account for non-perturbative effects. To compare to particle-level jet cross sections, supplementary 
calculations are required. PYTHIA and HERWIG-I— I- [l^] are used to generate samples without underlying 
event. Jets in these samples are reconstructed from partons after the parton shower, and observables are 
compared at the particle-level in the standard HERWIG-I--I- and PYTHIA samples. A multiplicative correction 
is calculated. This correction factor takes the NLO pQCD calculations to the particle level. The correction 
obtained using the PYTHIA AMBTl sample is taken as the default value for the analysis, and the systematic 
uncertainty is estimated from the maximum spread compared to the results from the other models. The size of 
this correction is less than 5% in all observables studied in the NLO pQCD analysis. The total uncertainty quoted 
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on the NLO pQCD calculations comes from the quadrature sum of the uncertainties from the renormalization 
and factorization scales, the proton pdfs, as and the non-perturbative corrections. 



5. Event Selection and Reconstruction 

A set of ATLAS first level (level-1) multi-jet triggers is used to select events for the analysis. Multi-jet triggers 
require several jets reconstructed with a level-1 sliding window algorithm. All multi-jet triggers are symmetric. 
Only two-jet and three-jet triggers were needed for the analysis. 

The primary vertex or vertices are found using tracks that originate in the beam collision spot 26], satisfy 
quality criteria [27] and have transverse momentum above 150 MeV. 

Topological calorimeter clusters evaluated at the electromagnetic scale pgj are used as inputs to the jet 
finding algorithm. The anti-fcf algorithm Q with resolution parameters R = 0.4 and R — 0.6 and full four- 
momentum recombination is used to reconstruct jets from clusters. The jet reconstruction is fully efficient in the 
Monte Carlo simulation for jets with transverse momentum above 30 GeV. The ATLAS Monte Carlo simulation 
compares well with the jet reconstruction efficiency measured with data p8| . 

Jets reconstructed at the electromagnetic scale are measured to have an energy which is lower than the true 
energy of interacting particles within the jet. This is caused by the different response of the calorimeter to 
hadrons and electrons and by energy lost in dead material, deposited outside the jet cone or excluded by the 
cell clustering algorithm. A Monte Carlo-based calibration that corrects for these effects as a function of 
and y is used to obtain jets with the correct energy scale [1^ . 

Jets considered in the analysis are selected using the following kinematic and data quality selection crite- 
ria: 

• An event must contain at least one jet with \y\ < 2.8 and a px greater than or equal to 80 GeV. 

• Jets are required to have |y| < 2.8 and px > 60 GeV in order to be counted. 

• A series of jet cleaning cuts were applied to eliminate various detector effects and suppress beam and other 
non-collision backgrounds. Overall, these cuts reduce the available statistics by less than 0.1% [301] . 

• Jets are only accepted if at least 70% of their charged particle pt comes from the event vertex. 

• Only events with a minimum of two selected jets are used in the analysis. 

For a total integrated luminosity of 2.43 pb"-'^, approximately 500,000 multi-jet events survived the selection 
cuts. 



6. Data Correction 

A correction is needed to compare the measurements to theoretical predictions. The correction, which accounts 
for trigger inefficiencies, detector resolutions and other detector effects that affect the jet counting, is performed 
in a single step using a bin-by-bin unfolding method calculated from Monte Carlo simulations. For each measured 
distribution, the corresponding Monte Carlo simulation cross section using truth jets as defined in Section |3] is 
evaluated in the relevant bins, along with the equivalent distributions obtained after the application of detector 
simulation and analysis cuts. The ratio of the true to the simulated distributions provides the multiplicative 
correction factor (unfolding factor) to be applied to the measured distributions. 

To perform the unfolding, ALPGEN-^HERWIG/ JIMMY AUETl Monte Carlo simulation is used. The un- 
folding uncertainty is estimated taking into account several effects. One arises from the spread in correction 
factors coming from different generators (ALPGEN+ HERWIG/JIMMY AUETl and PYTHIA AMBTl). A 
second detailed study is performed in which the simulated jet pr, y and 4> resolution is varied. Third, the shape 
of the simulated distributions is varied in order to account for possible biases caused by the input distributions. 
Samples with a trigger inefficiency in the crack region, with different pile-up rejection cuts and different primary 
vertex multiplicity distributions are also used to estimate the uncertainty arising from trigger effects and from 
the impact of overlapping proton-proton collisions. All these effects impact the unfolding systematics, and their 
uncertainties are ultimately added in quadrature to provide the final systematic uncertainty in the unfolding 
correction. Statistical uncertainties on the unfolding factors are important for certain bins and added to the 
total uncertainty. 

The combined systematic uncertainty on the data correction is less than 10% in most bins except those 
affected by high statistical uncertainties where the uncertainty can be up to 20%. 
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The systematic uncertainties in the luminosity calculation affect all cross section measurements, but cancel 
out in all measurements where cross-section ratios are involved. The luminosity of the dataset used in this 
paper has been calculated to be 2.43±0.08 pb~^ [Slj and the associated uncertainty is not shown in the figures. 

7. Uncertainty on the Jet Energy Scale 

The jet energy scale uncertainty is the dominant uncertainty for most results presented in this paper. The 
steeply falling cross sections as a function of jet implies that even a relatively small uncertainty in the 
determination of the jet translates into a substantial change in the cross sections as events migrate along 
the steeply falling curve. 

The jet energy scale in ATLAS [1^ has been calculated using jets from a dijet sample without near-by activity 
in the calorimeter. For a multi-jet analysis, a set of additional systematic uncertainties need to be considered. 
These uncertainties arise from changes in the calorimeter response to jets of different flavors as well as the 
presence of nearby activity in the calorimeter on the jet energy measurement. 

In summary, the jet energy scale uncertainty is primarily made of three components: the uncertainty calcu- 
lated for isolated jets, the uncertainty caused by the presence of nearby calorimeter deposits, and the flavor 
composition uncertainty. The uncertainty on the energy scale of isolated jets is the largest contributor to the 
total uncertainty in most bins, except for jets in the five and six-jet bins and of < 200GeV, for which the 
flavor composition uncertainty is comparable. The positive systematic uncertainty on the jet energy scale of 
isolated jets falling in the barrel and in high-multiplicity bins varies from 5% at 60 GeV to 2.5% at 1 TeV. In 
the three-jet and four-jet bins, where the flavor composition is better constrained, the systematic uncertainty is 
at most 3.5%. The negative systematic uncertainty is smaller and ^3% across all px in the barrel. The impact 
of nearby calorimeter deposits is small, increasing the overall uncertainty by at most 1%. 

8. Results 

In this section, measurements^ corrected to the particle level are compared to theoretical predictions. For 
comparisons to leading-order Monte Carlo simulations, the anti-fct algorithm with resolution parameter R = 0.4 
is used to define a jet. In Figures [T]-3 and 5(b), the orange error band bracketing the measured cross section 
corresponds to the total systematic uncertainty, evaluated by adding the individual systematic uncertainties in 
quadrature but excluding the uncertainty coming from the luminosity measurement. The ratio of the predictions 
from the Monte Carlo simulations to the measurements is shown at the bottom of each figure. All Monte Carlo 
simulations are normalized to the measured inclusive two jet cross section, these factors are shown on the 
relevant figures. 

Most ALPGEN Monte Carlo simulations studied predict an inclusive multi-jet cross section similar to the 
measured cross section, while the PYTHIA Monte Carlo simulation requires scaling factors which differ sig- 
nificantly from unity. The differences in the normalization factors between ALPGEN-I-PYTIIIA MC09' and 
ALPGEN-I-HERWIG/ JIMMY AUETl illustrate differences between PYTHIA and HERWIG/ JIMMY and their 
interplay with the matrix-element and parton-shower matching implemented in ALPGEN. The normalization 
factor for SHERPA is the closest to unity. 

Differences between the different ALPGEN tunes illustrate the large impact of the underlying event and parton 
shower tune on the leading-order predictions. It has been hypothesized that the different parton shower models 
available in PYTHIA and HERWIG may, through the MLM-matching used in ALPGEN, have a significant 
effect on the overall prediction for the cross sections. This issue would benefit from additional studies, and 
perhaps dedicated tuning of ALPGEN-I-PYTHIA or ALPGEN+HERWIG. 

Figure [IJa) shows the results for the cross section as a function of the inclusive jet multiplicity. The measure- 
ment systematics are dominated by the jet energy scale uncertainty and range from 10-20% at low multiplicities 
to almost 30-40% at high multiplicities. The Monte Carlo simulation predictions agree with the measured results 
across the full inclusive multiplicity spectrum. 

A study that reduces significantly the impact of systematic uncertainties is the ratio of the n-jet to (n — 1)- 
jet cross section as a function of multiplicity. In this ratio, the impact of the jet energy scale uncertainty is 



■^AU measurements in this section have been compiled in Tables that can be found in HEPDATA. The NLO pQCD calculation 
results are also presented in the Tables when applicable. 
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significantly reduced and the uncertainty due to the luminosity cancels out. Figure [Ijb) shows the results for 
such a study. 

Both the data correction and the jet energy scale uncertainties contribute comparably to the total systematic 
uncertainty, whereas the statistical uncertainties are smaller than the systematic uncertainties, and negligible 
in most bins. All Monte Carlo simulations agree well with the data, yet there is a noticeable spread of ^-^15-20% 
in their predictions. 

The differential cross section for multi-jet events as a function of the jet is useful for characterizing the 
kinematic features. Figure [2] presents the px-dependent differential cross sections for the leading and third 
leading jet in multi-jet events. The systematic uncertainty in the measurement is 10-20% across and increas- 
ing up to 30% for the fourth leading jet differential cross section. The jet energy scale systematic uncertainty 
remains the dominant uncertainty in the measurement. All Monte Carlo simulations agree reasonably well with 
the data. The PYTHIA AMBTl Monte Carlo simulation predicts a somewhat steeper slope, compared to the 
data, as a function of the leading jet and of the second leading jet px- The effect is most noticeable in the 
highest px bins. 
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Figure 1: Total inclusive jet cross section as a function of multiplicity (a). Ratio of the n-jet cross section to the (n— 1) 
cross section for values of n varying from three to six (b). 
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Figure 2: Differential cross section as a function of leading jet pt for events with A'jcts > 2 (a), S'^'' leading jet pT for 
events with A'jcts > 3 (b). 
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Figure 3: Three-to-two-jet differential cross-section ratio as a function of tlie leading jet pt- In tlie figures, a resolution 
parameter R = 0.6 is used. Tlie three figures contain a minimum pt cut for all non-leading jets of (a) 60 GeV (b) 80 GeV 
and (c) 110 GeV. 
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Figure 4: Three-to- two-jet differential cross-section ratio as a function of the leading jet px- In the figures a resolution 
parameter R — 0.6 is used. The three figures contain a minimum px cut for all non-leading jets of (a) 60 GeV (b) 
80 GeV and (c) 110 GeV. The resuhs are compared to a NLO pQCD calculation with the MSTW 2008 nio pdf set . 
The systematic uncertainties on the theoretical prediction are shown as dotted red lines above and below the theoretical 
prediction. 



A measurement with particular sensitivity to limitations in the leading-order Monte Carlo simulations and 
NLO pQCD calculations is the ratio of the inclusive three-to- t-wo-jet differential cross section as a function of 
some characteristic scale in the event. In this measurement, the uncertainty in the luminosity determination 
cancels out, uncertainties in the jet energy scale are reduced, and statistical uncertainties are limited only by 
the inclusive three-jet sample. 

For tuning Monte Carlo simulations, the three-to-two-jet ratio as a function of the leading jet has been 
shown to be particularly sensitive to effects from final state radiation [s^] ■ Figure [3] presents the results on the 
measurement of the threc-to-two-jet cross section ratio as a function of leading jet for jets built with the 
anti-fcf algorithm using the resolution parameter R=0.6 and with different minimum cuts for all non-leading 
jets^. The systematic uncertainties on the measurement are small (^5%), except in the lowest px bin, where 
uncertainties in the unfolding correction dominate. ALPGEN+HERWIG AUETl and ALPGEN-I-PYTHIA 
MC09' describe the data well, and the agreements appear to be largely independent of the tunes chosen. 
SHERPA also describes the data well. PYTHIA AMBTl predicts a higher ratio than what is measured over the 
Pt range from 200 GeV to 600 GeV^. The systematic uncertainty in the lowest pr bin decreases significantly 
as the minimum px cut is raised to 80 GeV for all jets. 



^Results (not shown) were also obtained using R = 0.4 and are compiled in Tables in HEPDATA. 

^Further results for more Monte Carlo simulations are compiled in Tables in HEPDATA. The disagreement observed when 
PYTHIA Perugia2010, PYTHIA MC09 and HERWIG-|-|- are used is similar to that observed here. 
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Figure 5: Three-to-two-jet differential cross-section ratio as a function of tlie sum of the pt of the two leading jets using 
R = 0.6. The left and right figures present the same measurements and error bands. The data are compared to (a) 
a NLO pQCD calculation and (b) several leading-order Monte Carlo simulations. The systematic uncertainties on the 
theoretical prediction for the NLO pQCD calculations are shown as dotted red lines above and below the theoretical 
prediction. 



Figure 3] presents the same measurement results as Figure [31 except the data are now compared to the NLO 
pQCD calculations corrected for non-perturbative effects. The MSTW 2008 nlo pdf set has been used, but 
comparable results are obtained with the CTEQ 6.6 pdf set. The systematic uncertainties on the theoretical 
predictions are shown as dotted red lines above and below the theoretical prediction. The NLO pQCD calcu- 
lations describe the data well, except in the lowest bin, where there is a large discrepancy. The discrepancy 
diminishes significantly once the minimum px for all jets is raised to 110 GeV and the px of the leading jet is 
required to be greater than 160 GeV. 

Additional NLO pQCD calculations of the three-to-two-jet cross section ratio were performed as a function of 
different kinematic variables, such as Ht , the sum of the px of the two leading jets {Hj. ) and the sum of the px 

(2) 

of the three leading jets. The NLO pQCD calculation for the ratio as a function of was found to give the 
smallest theoretical scale uncertainty and is, therefore, most sensitive to input parameters such as as- Figure [5| 
shows a comparison of the measurement to both NLO pQCD (a) and leading order calculations (b) for R = 
0.6. Scale uncertainties of the NLO pQGD calculations are larger for jets with R = 0.4 than with R = 0.6. The 
theoretical uncertainty of the NLO pQGD calculations shown in Figure |4| is comparable to the measurement 
uncertainties, but is significantly reduced compared to theoretical uncertainties presented in Figure [S] With the 
reduced theoretical uncertainty, the disagreement between data and the NLO pQCD calculations in the lowest 

(2) 

Hjp bin is now enhanced. Due to the kinematic cuts applied in the analysis, the NLO pQCD calculations only 
account for the lowest-order contribution to the two-jet cross section in the region where the sum of the first and 
second leading jet px is less than 160 GeV. Consequentially, this effective leading-order estimation is subject to 
large theoretical uncertainties, which might be responsible for the observed discrepancy. 

A comparison of the same measurement to leading-order Monte Carlo simulations is given in Figure [Sj (b). 
The general agreement between leading-order Monte Carlo simulations with the measurements follows the same 
general trends as the comparison of the three-to- two-jet ratio versus leading jet px shown in Figure [D 



9. Summary and Conclusion 

A first dedicated study of multi-jet events has been performed using the ATLAS detector at a center-of-mass 
energy of 7 TeV with an integrated luminosity of 2.43 pb~^. For events containing two or more jets with 
px > 60 GeV, of which at least one has px > 80 GeV, good agreement is found between data and leading-order 
Monte Carlo simulations with parton-shower tunes that describe adequately the ATLAS yfs = 7 TeV underlying- 
event data. The agreement is found after the predictions of the Monte Carlo simulations are normalized to the 
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measured inclusive two-jet cross section. The models have been compared to multi-jet inclusive and differential 
cross sections. The present study extends up to a multiplicity of six jets, up to jet px of 800 GeV and up to 
event of 1.6 TeV. 

All models reproduce the main features of the multijet data. ALPGEN, which contains higher multiplicity 
tree-level matrix elements matched to parton showers, generally describes the shapes well, whether the HERWIG 
or PYTHIA parton showers are used. Similarly, SHERPA describes the data well. However, PYTHIA, which 
contains a 2 — )■ 2 leading-order matrix element augmented by parton showers, shows only a marginal agreement 
with the data. This is most noticeable in the three-to- two-jet cross section ratios. 

A measurement of the three-to- two-jet cross section ratio as a function of the leading jet and the sum of 
the two leading jet pxs has also been performed and is described by ALPGEN, SHERPA and a NLO pQCD 
calculation, albeit with a significant discrepancy in the lowest bin for the latter comparison. Comparisons 
with NLO pQCD calculations may be useful for constraining parameters, such as parton distribution functions 
or the value of the strong coupling constant, as, given that the systematic uncertainties from the measurement 
are comparable to the theoretical uncertainties. 
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